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Abstract 

We calculate the partition function Z(t) and the asymptotic integrated level density N(E) for 
Yang-Mills-Higgs Quantum Mechanics for two and three dimensions (n = 2, 3). Due to the infinite 
volume of the phase space T on energy shell for n = 2, it is not possible to disentangle completely 
the coupled oscillators (x 2 y 2 -model) from the Higgs sector. The situation is different for n = 3 
for which T is finite. The transition from order to chaos in these systems is expressed by the 
corresponding transitions in Z{t) and N{E), analogous to the transitions in adjacent level spacing 
distribution from Poisson distribution to Wigner-Dyson distribution. We also discuss a related 
system with quartic coupled oscillators and two dimensional quartic free oscillators for which, 
contrary to YMHQM, both coupling constants are dimensionless. 

PACS numbers: 05.45.Mt, ll.15.-q, 03.65.Sq, 11.15.Kc 



* Present address: 3106 Hornbuckle Place, Durham, NC 27707, USA. E-mail: sgmatin@aol.com 

• E-mail: yjng@physics.unc.edu 



1 



I. INTRODUCTION 

The discovery of chaoticity in the classical Yang-Mills (YM) equations [j| has attracted 
attention to the system of coupled quartic oscillators with the potential x 2 y 2 , where x and y 
are functions of time only This system is the simplest limiting case for the homogeneous YM 
equations (the so-called YM Classical Mechanics) with n = 2 degrees of freedom. Despite its 
simplicity, the x 2 y 2 model exhibits a rich versatile chaotic behaviour and belongs to the most 
chaotic potential systems known. Not surprisingly, this potential has been used in many 
fields, including chemistry, astronomy, astrophysics, and cosmology (chaotic inflation). 

From the quantum mechanical point of view, this system (YM Quantum Mechanics) has 
a discrete spectrum 0, [J in spite of having an infinite volume of energetically accessible 



phase space T 



r = J dxdydxdy ^i^x 2 + i^J 2 + i^ x2 y 2 ~~ E 



where the dot stands for d/dt. This potential, therefore, violates the Weil's law[y|, a semi- 
classical relation, which states that the average number N(E) of quantum energy levels with 
energy less than E is asymptotically proportional to T. 

Some time ago, there appeared two important papers[?], 0] devoted to the calculation 
of the partition function Z(t) and the asymptotic integrated level density N(E) for the 
x 2 y 2 model. Their method was based on an adiabatic separation in the partition function's 
dependence on x and y out in the narrow channels of the equipotential surface \xy\ = 
constant. Remarkably the dependence on the boundary dividing the two regions (central 
(\x\ < Q, \y\ < Q) and channel (Q < \x\, \y\ < oo)) with quite different physics (essentially 
classical for the first region and intrinsically quantum for the second region) disappears 
in the final result for the partition function. This insensitivity to the dependence on the 
boundary (the value of Q) may seem to bode well for this method of calculation. 

However, a priori we do not expect that the partition function of a non-integrable system 
with infinite phase space volume is calculable without any approximation. Therefore we seek 
an alternative approach which brings to fore the fact that T = oo for the n = 2 case. With 
this in mind, in this paper we calculate the partition function and the asymptotic integrated 
level density for the so-called Yang-Mills-Higgs Quantum Mechanics (YMHQM) jj, |J It is 
interesting to examine how Z(t) and N(E) behave in the limit the Higgs coupling to the x— 
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and y— amplitudes vanishes and one recovers the pure x 2 ?/ 2 -system. We also consider the 
YMHQM for the case of n = 3 for which T is finite 0, 0] even for the pure YM system. 

II. YANG-MILLS-HIGGS MECHANICS 

There are several mechanisms that can suppress the classical chaoticityof chaotic systems 
with dimension n = 2 and, in particular, to the the YM system (see Ref. 5|). One of them is 
the well-known Higgs mechanism. For spatially homogeneous fields, if only the interaction 
of the YM fields with the Higgs vacuum is considered, the classical Hamiltonian density for 
n = 2 is given by 

H= 1 -(^ + f) + 9 ^xV + v ^ 2 + y 2 )- (i) 

It is known that there is a classical "phase transition" from chaos to regular motions 
as the vacuum expectation value of the Higgs field (0) = v gets large enough. Here there 

2 4 

is one control parameter k = a -g-, where H is the conserved energy density. At large k the 
motion is regular, whereas at k = one deals with the developed chaos of YM Classical 
Mechanics. In fact, chaos appears already at k ~ 0.60. [9] Therefore, in the study of the 
quantum counterpart of Eq. (^), we expect that there is a transition from one type of N(E) 
to another, depending on the parameter v. The analogous transition in the adjacent energy 
level s pac ing distribution, as a function of v, was predicted jl^] and established in several 
papers 

In the next section we calculate the partition function for the quantum-mechanical coun- 
terpart of the classical Hamiltonian density Eq. In section IV, we calculate the in- 
tegrated level density. Section V is devoted to the n = 3 case. We conclude in the 
last section with some discussions involving a related system with potential of the form 

^V + ?(* 4 + z/ 4 )- 

III. PARTITION FUNCTION FOR YMH QUANTUM MECHANICS 

The quantum Hamiltonian corresponding to Eq. (Q) is given by (mass m = 1 in the 
following) 

H = ^ + 9 4x 2 y 2 + V 4^ + y*), (2) 
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where is the momentum operator. First a word on units. All quantities below 

are given in units of energy E: [H] = E, [t] = E~\ [x] = [y] = E 1 / 4 , [g] = E°, [v] = E 1 / 4 , 
[h] = E 3 ^ 4 . It is obvious that the operator Eq. (J2J) has a discrete spectrum. 

Here we calculate the trace of the heat kernel exp(— tH), the partition function for the 
Hamiltonian operator. The trace is defined for any quantum operator Aw in the Wigner 
representation (i.e., in the classical 2n- dimensional phase space) 

Tr(Aw) = ( 27r fi)n / d'pd'^AwC^, qf). (3) 

Thus, we have 

Z{t) = Tr{e~ tn ) = / dEe~ tE p(E). (4) 
Jo 

The second equation expresses the fact that the partition function Z(t) and the density of 
eigenstates p(E) form a Laplace transform pair. The asymptotic integrated density of states 
N(E) is given by the inverse Laplace transform of Z(t)/t 

N(E)= f E dE'p(E') = L-'(^P]. (5) 



t 

Integrating Eq. (JIJ over p x and p y for the Hamiltonian given by Eq. (J2J), we get 

Z (t) = t-^j— [ dxdye-^^^ 2 ^ 2 . (6) 



(2nh) 2 t 

After the integration over x, we can make a change of variable from y to s — ^tv 2 y 2 to 
obtain 

Z(t) = ^ (^)"\^ fte-.-V> fl + A)-" 2 , (7) 
v ; (2nh) 2 \ t ) \ tv 4 Jo V 2z) K ' 



with z = ^2 ■ But aside from a factor of y/2ze z , the integral is just a representation of ifo(z), 
the modified Bessel function of the third kind or the MacDonald function of order 0. (see, 
e.g., p. 140 of Ref.Q]). Thus we arrive at the precise expression for the partition function 
corresponding to the Hamiltonian Eq. (J2J) 



1 1 (tv 4 \ T ^ (tv 4 \ 



So far we have not made any calculational approximations. For two uncoupled oscillators 
(corresponding to g — 0) Eq. (jHJ) gives (with the aid of K (z) ~ \nt e ~ z as z "~ > °°) 



as it should. For the sake of completeness, let us also give the partition function for one 
oscillator (n — 1) 

From Eq. (|SJ), it is obvious that Z(t) diverges logarithmically in the limit v — > (the 
x 2 y 2 model). Indeed, using Kq(z) ~ — log(z/2) — C (with C being the Euler constant) as 
z — ► 0, we get, for t> ~ 0, 

z(t ^vb^( log SH' (10) 

The impossibility to disentangle the coupled oscillators from the uncoupled ones is a reflec- 
tion of the logarithmic divergence of the phase space volume V on the .E-shell at v — for 
n = 2, as we have mentioned above. For n = 3 (see section V below), T is finite and the 
precise v — limit exists. 



IV. INTEGRATED LEVEL DENSITY iV(.E) (ra = 2) 

It is well-known that the asymptotic integrated density of states N(E) is related to the 
small t divergence of the partition function Z(t), according to the Karamata-Tauberian 
theorem (see, e.g., Ref. 3]). Using Eq. (jHJ) and Eq. (jSJ), we obtain the precise expression 
for N(E) (with the aid of formula 3.16.3.3 of Ref.pj) 



4 ' 



E 2 (\ 1 2g 2 E' 



~ 2nV F V2' 2 ;3; ^ J' 

where F is the Gauss hypergeometric function. For g = (the case of two uncoupled 
oscillators) we have, from Eq. (|TT| . 

N(E) = ^~ 2 E\ (12) 
in agreement with Eq. (jH)). For completeness, for one oscillator, we have, from Eq. (pFj) . 

jV(E) = — E. (12') 

AW 
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For the v — > limit, one needs the asymptotic expression of F(a, b; c; z). Using formula 



9.7.7 of Ref. 



12J, we get 

2^/2 E^ 2 /. g 2 E 

^WK 0g ^~-— 3 



N(E) ~ -i^-^r ( log y — + 5 log 2 - - ] . (13) 



One can verify that the Z(t) and N(E) for Eq. (JJJJ), Eq. (|FJ), and the corresponding Eq. (fT2j) . 
Eq. ([12'p . and the ratio of the logarithmic parts of Eq. ()10|) and Eq. (|13|). are in accordance 
with Theorems 1.1 and 1.4 of Ref. j^. Also from Eq. ()12j) and Eq. ()13|) . we see that the 
transition from order (g = 0) to chaos (g ^ and large, v small, corresponds to the 

change in the E'-dependence of N(E) from N(E) ~ E 2 to N(E) ~ i? 3 / 2 log £7, analogous 
to the change in the neighbour level spacing distribution from the Poisson distribution to 
the Wigner-Dyson distribution. It can be seen that as the power a of the homogeneous 
potential \xy\ a increases from 1 to oo, there is [3] a systematic decrease in the power of 
E in N(E) from 2 to 1, the latter corresponding to the case of the hyperbola billiard with 
N(E) ~ E logE (see Ref. 3|). This observation may be relevant to the Hilbert-Polya-Berry 
program in identifying a quantum (chaotic) Hamiltonian whose eigenvalues reproduce the 
Riemann zeta-function zeros. 



V. SOME REMARKS ON THE THREE DIMENSIONAL YM AND YMH QUAN- 
TUM MECHANICS 

Let us now consider the YMH Quantum Mechanics for n = 3 with the Hamiltonian (with 
mass m = 1 again) 

n, L — >2 i 1/ / 2 2 i 22, 2 2\ , " / 2 2, 2\ /i a\ 

H = -p +^(xy +y z +zx) + —(x +y +z). (14) 

For the partition function, the integrations over p? can be easily carried out. Introducing 
the cylindrical coordinates (x = r cos <f),y = r sin <f>,z = z) and performing the integrations 
over z and 0, we obtain 

1 /2tt\ 2 f°° rdr ( tv 2 r 2 tg 2 r A \ (tg 2 r A \ . . 

z[t) = (2^p It J 271 L ^wrw exp l"^ ^) °V^) ' (15) 

where Iq(z) is the modified Bessel function of the first kind of order 0. 



1 We need to emphasize that, from the results of Ref.|9(, the chaotic regime for N(E) actually begins at 
finite v = {nE/g 2 ) 1 / 4 . 
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We see that the integral in Eq. (JT5j) is well behaved at v = 0; it is equal to 



The integration over u can be done by using formula 2.15.3.3 of Ref.[3] to yield 

r 3 (|) i 

Z ^ = 2 7/4 7] 4 3/2 gWtftflV ( 16 ) 

which agrees with Eq. (4.5) of Ref.|8[ (after taking into account the extra factor of | in the 
definition of the potential in Eq. (|14|)). The difference between the n = 2 and n = 3 cases is 
due to the fact that for n = 3 the energetically accessible region pinches as 1/x 2 at large x, 
in contrast to the n = 2 case for which it pinches as 1/x, leading to the logarithmic divergent 
result (Eq. ([10)1 ) at v — 0. Note that the contibution from the channels {t l ^x » 1) for 
n = 3 is negligible in the method of separating the central part from the channels used in 
Ref. 0, Isj]. Thus for n = 3 the coupled oscillators are disentangled from the free oscillators. 
From Eq. (|lfij). it follows that 

which agrees with the result given in Ref.Q (the first of Eq. (5.19)). 

Finally, for the sake of completeness, here are the expressions for Z(t) and N(E) at g = 
for the case of n = 3. Eq. (fTHjl yields 

Z(() = (^F- (18) 

Of course, Eq. 0, Eq. (jO 7 )) . and Eq. (fTKj) are in accordance with the well-known relation 
of the trace of the n-harmonic-oscillator heat kernel, viz., Z(t) = (2sinh^) in the 

hn v 2 J 

quasiclassical approximation (t 0)[7j,l3j. From Eq. (|T5j) we have, for N(E), 

N ^ = T^E 3 - (19) 

DAI f 1:1 

Eq. (|Td^ - Eq. (|19|) show that the Z(t) and N(E) for the n = 3 case are in full agreement 
with Theorem 1.1 of Ref. [3j]. 



VI. DISCUSSIONS 



From the previous two sections, we see that the log-E factor in N(E) is specific to the 
x 2 y 2 -mode\ (with or without the Higgs term), and is not generic to chaotic systems with 
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potentials \xy\ a (a = 2,3,...) for n > 2. But one may wonder whether it is generic to 
other homogeneous potentials for n = 2. The answer depends on the type of homogeneous 
potentials added to the x 2 y 2 term. Consider the following example of n = 2 system (with 
finite phase space), described by the following quantum Hamiltonian 

W=^ + y*V+ j(x 4 + y 4 ). (20) 

Like g, b is dimensionless. It is known that for b 2 = g 2 /3 the system Eq. (|2T)|) is integrable. 
In the calculation of Z(t), the integrations over x and y can best be done by using the polar 
coordinates. After the gaussian integration over the radial coordinate, the integration over 



the angle yields the complete elliptic integral of the first kind 
Gauss hypergeometric function. The result is given by 



121 ] which is related to the 



V^r 1 (I I b 2 -g 2 \ 
The corresponding asymptotic integrated density of states N(E) is given by 



For the b — > limit, we use formula 9.7.7 of Ref . |12j again to get 



z(t) -7SW^ los ^ (23) 



We note the logarithmic divergence as b — > 0, i.e., when one tries to disentangle the coupled 
quartic oscillators from the free quartic oscillators. A simple inverse Laplace transform yields 
the level density 

Just as there is no logt dependence in Z(t), we find no \ogE dependence in N(E) for this 
system. This result is not surprising from the viewpoint of dimensional considerations since 
both couplings in this example are dimensionless. In the case of YMHQM, the coupling v is 
dimensionful, and due to this, log t and log E appear in Eq. (fTUjl and Eq. (fTHjl respectively 
in the v — > limit. 

For completeness we also give the partition function for the g = case which corresponds 
to an integrable system. Using F(~,|;l;|) = r 2 (|)/(27r 3 / 2 ), we have 

Z(t) = UJ 9 - . (25) 
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The corresponding asymptotic integrated density of states is given by 

r 2 fl) £3/2 

For another integrable case with b 2 = g 2 /3 in this examplej3], using ~; 1; — 1) = 

r 2 (i)/(2^) 3 / 2 , we get 

z «> = ^spW- (27) 

and 

Y 2 (l) £3/2 

ME) = , \i /0 9 . (28) 
v ; 3(2tt 3 ) 1 /2 bh 2 V ; 

Finally, for the third integrable case with b = g we immediately obtain the Z(t) and N(E) 

from Eq. (0} and Eq. (J22J) using |; 1,0) = 1. 
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